Marc van Regenmortel, personal recollections on a forward-thinking editor.

Marc van Regenmortel was the Editor-in-Chief of the Journal of Molecular Recognition for the last 25 years. Without attempting to summarize Marc's exceptional career and achievements, we would like to tell the story of the tortuous and contingent path to the unravelling of a key molecular recognition process in antigenicity. Life is indeed full of contingencies and scientific life, full of meetings and random encounters, is prone to contingencies, a key element in discovery and innovation.

Contribution of tRNA sequence and modifications to the decoding preferences of E. coli and M. mycoides tRNAGlyUCC for synonymous glycine codons.

tRNA superwobbling, used by certain bacteria and organelles, is an intriguing decoding concept in which a single tRNA isoacceptor is used to decode all synonymous codons of a four-fold degenerate codon box. While Escherichia coli relies on three tRNAGly isoacceptors to decode the four glycine codons (GGN), Mycoplasma mycoides requires only a single tRNAGly. Both organisms express tRNAGly with the anticodon UCC, which are remarkably similar in sequence but different in their decoding ability. By systematically introducing mutations and altering the number and type of tRNA modifications using chemically synthesized tRNAs, we elucidated the contribution of individual nucleotides and chemical groups to decoding by the E. coli and M. mycoides tRNAGly. The tRNA sequence was identified as the key factor for superwobbling, revealing the T-arm sequence as a novel pivotal element. In addition, the presence of tRNA modifications, although not essential for providing superwobbling, was shown to delicately fine-tune and balance the decoding of synonymous codons. This emphasizes that the tRNA sequence and its modifications together form an intricate system of high complexity that is indispensable for accurate and efficient decoding.

SilkMeta: a comprehensive platform for sharing and exploiting pan-genomic and multi-omic silkworm data.

The silkworm Bombyx mori is a domesticated insect that serves as an animal model for research and agriculture. The silkworm super-pan-genome dataset, which we published last year, is a unique resource for the study of global genomic diversity and phenotype-genotype association. Here we present SilkMeta (http://silkmeta.org.cn), a comprehensive database covering the available silkworm pan-genome and multi-omics data. The database contains 1082 short-read genomes, 546 long-read assembled genomes, 1168 transcriptomes, 294 phenotype characterizations (phenome), tens of millions of variations (variome), 7253 long non-coding RNAs (lncRNAs), 18 717 full length transcripts and a set of population statistics. We have compiled publications on functional genomics research and genetic stock deciphering (mutant map). A range of bioinformatics tools is also provided for data visualization and retrieval. The large batch of omics data and tools were integrated in twelve functional modules that provide useful strategies and data for comparative and functional genomics research. The interactive bioinformatics platform SilkMeta will benefit not only the silkworm but also the insect biology communities.

Assessment of three-dimensional RNA structure prediction in CASP15.

The prediction of RNA three-dimensional structures remains an unsolved problem. Here, we report assessments of RNA structure predictions in CASP15, the first CASP exercise that involved RNA structure modeling. Forty-two predictor groups submitted models for at least one of twelve RNA-containing targets. These models were evaluated by the RNA-Puzzles organizers and, separately, by a CASP-recruited team using metrics (GDT, lDDT) and approaches (Z-score rankings) initially developed for assessment of proteins and generalized here for RNA assessment. The two assessments independently ranked the same predictor groups as first (AIchemy_RNA2), second (Chen), and third (RNAPolis and GeneSilico, tied); predictions from deep learning approaches were significantly worse than these top ranked groups, which did not use deep learning. Further analyses based on direct comparison of predicted models to cryogenic electron microscopy (cryo-EM) maps and x-ray diffraction data support these rankings. With the exception of two RNA-protein complexes, models submitted by CASP15 groups correctly predicted the global fold of the RNA targets. Comparisons of CASP15 submissions to designed RNA nanostructures as well as molecular replacement trials highlight the potential utility of current RNA modeling approaches for RNA nanotechnology and structural biology, respectively. Nevertheless, challenges remain in modeling fine details such as noncanonical pairs, in ranking among submitted models, and in prediction of multiple structures resolved by cryo-EM or crystallography.

SARS-CoV-2 NSP1 induces mRNA cleavages on the ribosome.

In severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the non-structural protein NSP1 inhibits translation of host mRNAs by binding to the mRNA entry channel of the ribosome and, together with the 5'-untranslated region (UTR) of the viral mRNAs, allows the evasion of that inhibition. Here, we show that NSP1 mediates endonucleolytic cleavages of both host and viral mRNAs in the 5'UTR, but with different cleavage patterns. The first pattern is observed in host mRNAs with cleavages interspersed regularly and close to the 5' cap (6-11 nt downstream of the cap). Those cleavage positions depend more on the position relative to the 5' cap than on the sequence itself. The second cleavage pattern occurs at high NSP1 concentrations and only in SARS-CoV-2 RNAs, with the cleavages clustered at positions 45, 46 and 49. Both patterns of cleavage occur with the mRNA and NSP1 bound to the ribosome, with the SL1 hairpin at the 5' end sufficient to protect from NSP1-mediated degradation at low NSP1 concentrations. We show further that the N-terminal domain of NSP1 is necessary and sufficient for efficient cleavage. We suggest that in the ribosome-bound NSP1 protein the catalytic residues of the N-terminal domain are unmasked by the remodelling of the α1- and α2-helices of the C-terminal domain.

Assessment of three-dimensional RNA structure prediction in CASP15.

The prediction of RNA three-dimensional structures remains an unsolved problem. Here, we report assessments of RNA structure predictions in CASP15, the first CASP exercise that involved RNA structure modeling. Forty two predictor groups submitted models for at least one of twelve RNA-containing targets. These models were evaluated by the RNA-Puzzles organizers and, separately, by a CASP-recruited team using metrics (GDT, lDDT) and approaches (Z-score rankings) initially developed for assessment of proteins and generalized here for RNA assessment. The two assessments independently ranked the same predictor groups as first (AIchemy_RNA2), second (Chen), and third (RNAPolis and GeneSilico, tied); predictions from deep learning approaches were significantly worse than these top ranked groups, which did not use deep learning. Further analyses based on direct comparison of predicted models to cryogenic electron microscopy (cryo-EM) maps and X-ray diffraction data support these rankings. With the exception of two RNA-protein complexes, models submitted by CASP15 groups correctly predicted the global fold of the RNA targets. Comparisons of CASP15 submissions to designed RNA nanostructures as well as molecular replacement trials highlight the potential utility of current RNA modeling approaches for RNA nanotechnology and structural biology, respectively. Nevertheless, challenges remain in modeling fine details such as non-canonical pairs, in ranking among submitted models, and in prediction of multiple structures resolved by cryo-EM or crystallography.

Mistranslation of the genetic code by a new family of bacterial transfer RNAs.

The correct coupling of amino acids with transfer RNAs (tRNAs) is vital for translating genetic information into functional proteins. Errors during this process lead to mistranslation, where a codon is translated using the wrong amino acid. While unregulated and prolonged mistranslation is often toxic, growing evidence suggests that organisms, from bacteria to humans, can induce and use mistranslation as a mechanism to overcome unfavorable environmental conditions. Most known cases of mistranslation are caused by translation factors with poor substrate specificity or when substrate discrimination is sensitive to molecular changes such as mutations or posttranslational modifications. Here we report two novel families of tRNAs, encoded by bacteria from the Streptomyces and Kitasatospora genera, that adopted dual identities by integrating the anticodons AUU (for Asn) or AGU (for Thr) into the structure of a distinct proline tRNA. These tRNAs are typically encoded next to a full-length or truncated version of a distinct isoform of bacterial-type prolyl-tRNA synthetase. Using two protein reporters, we showed that these tRNAs translate asparagine and threonine codons with proline. Moreover, when expressed in Escherichia coli, the tRNAs cause varying growth defects due to global Asn-to-Pro and Thr-to-Pro mutations. Yet, proteome-wide substitutions of Asn with Pro induced by tRNA expression increased cell tolerance to the antibiotic carbenicillin, indicating that Pro mistranslation can be beneficial under certain conditions. Collectively, our results significantly expand the catalog of organisms known to possess dedicated mistranslation machinery and support the concept that mistranslation is a mechanism for cellular resiliency against environmental stress.

Anionic G•U pairs in bacterial ribosomal rRNAs.

Wobble GU pairs (or G•U) occur frequently within double-stranded RNA helices interspersed between standard G=C and A-U Watson-Crick pairs. Another type of G•U pair interacting via their Watson-Crick edges has been observed in the A site of ribosome structures between a modified U34 in the tRNA anticodon triplet and G + 3 in the mRNA. In such pairs, the electronic structure of the U is changed with a negative charge on N3(U), resulting in two H-bonds between N1(G)…O4(U) and N2(G)…N3(U). Here, we report that such pairs occur in other highly conserved positions in ribosomal RNAs of bacteria in the absence of U modification. An anionic cis Watson-Crick G•G pair is also observed and well conserved in the small subunit. These pairs are observed in tightly folded regions.

Identification of Up47 in three thermophilic archaea, one mesophilic archaeon, and one hyperthermophilic bacterium.

Analysis of the profile of the tRNA modifications in several Archaea allowed us to observe a novel modified uridine in the V-loop of several tRNAs from two species: Pyrococcus furiosus and Sulfolobus acidocaldarius Recently, Ohira and colleagues characterized 2'-phosphouridine (Up) at position 47 in tRNAs of thermophilic Sulfurisphaera tokodaii, as well as in several other archaea and thermophilic bacteria. From the presence of the gene arkI corresponding to the RNA kinase responsible for Up47 formation, they also concluded that Up47 should be present in tRNAs of other thermophilic Archaea Reanalysis of our earlier data confirms that the unidentified residue in tRNAs of both P. furiosus and S. acidocaldarius is indeed 2'-phosphouridine followed by m5C48. Moreover, we find this modification in several tRNAs of other Archaea and of the hyperthermophilic bacterium Aquifex aeolicus.

Data, data, burning deep, in the forests of the net.

Continuous and imaginative technological developments are leading to a massive accumulation of various types of data in all areas of biological research. As a result, the central importance of databases is increasing. Databases related to biology must not only be structured using controlled vocabularies, but also be fully integrated into the whole biological domain. To achieve this goal, they must be systematically grounded in biological evolution and exploit the available tools of evolutionary systematics to contribute to our understanding of life processes.

Computational Pipeline for Reference-Free Comparative Analysis of RNA 3D Structures Applied to SARS-CoV-2 UTR Models.

RNA is a unique biomolecule that is involved in a variety of fundamental biological functions, all of which depend solely on its structure and dynamics. Since the experimental determination of crystal RNA structures is laborious, computational 3D structure prediction methods are experiencing an ongoing and thriving development. Such methods can lead to many models; thus, it is necessary to build comparisons and extract common structural motifs for further medical or biological studies. Here, we introduce a computational pipeline dedicated to reference-free high-throughput comparative analysis of 3D RNA structures. We show its application in the RNA-Puzzles challenge, in which five participating groups attempted to predict the three-dimensional structures of 5'- and 3'-untranslated regions (UTRs) of the SARS-CoV-2 genome. We report the results of this puzzle and discuss the structural motifs obtained from the analysis. All simulated models and tools incorporated into the pipeline are open to scientific and academic use.

Eukaryotic tRNA sequences present conserved and amino acid-specific structural signatures.

Metazoan organisms have many tRNA genes responsible for decoding amino acids. The set of all tRNA genes can be grouped in sets of common amino acids and isoacceptor tRNAs that are aminoacylated by corresponding aminoacyl-tRNA synthetases. Analysis of tRNA alignments shows that, despite the high number of tRNA genes, specific tRNA sequence motifs are highly conserved across multicellular eukaryotes. The conservation often extends throughout the isoacceptors and isodecoders with, in some cases, two sets of conserved isodecoders. This study is focused on non-Watson-Crick base pairs in the helical stems, especially GoU pairs. Each of the four helical stems may contain one or more conserved GoU pairs. Some are amino acid specific and could represent identity elements for the cognate aminoacyl tRNA synthetases. Other GoU pairs are found in more than a single amino acid and could be critical for native folding of the tRNAs. Interestingly, some GoU pairs are anticodon-specific, and others are found in phylogenetically-specific clades. Although the distribution of conservation likely reflects a balance between accommodating isotype-specific functions as well as those shared by all tRNAs essential for ribosomal translation, such conservations may indicate the existence of specialized tRNAs for specific translation targets, cellular conditions, or alternative functions.

Correlated sequence signatures are present within the genomic 5'UTR RNA and NSP1 protein in coronaviruses.

The 5'UTR part of coronavirus genomes plays key roles in the viral replication cycle and translation of viral mRNAs. The first 75-80 nt, also called the leader sequence, are identical for genomic mRNA and subgenomic mRNAs. Recently, it was shown that cooperative actions of a 5'UTR segment and the nonstructural protein NSP1 are essential for both the inhibition of host mRNAs and for specific translation of viral mRNAs. Here, sequence analyses of both the 5'UTR RNA segment and the NSP1 protein have been done for several coronaviruses, with special attention to the betacoronaviruses. The conclusions are: (i) precise specific molecular signatures can be found in both the RNA and the NSP1 protein; (ii) both types of signatures correlate between each other. Indeed, definite sequence motifs in the RNA correlate with sequence motifs in the protein, indicating a coevolution between the 5'UTR and NSP1 in betacoronaviruses. Experimental mutational data on 5'UTR and NSP1 from SARS-CoV-2 using cell-free translation extracts support these conclusions and show that some conserved key residues in the amino-terminal half of the NSP1 protein are essential for evasion to the inhibitory effect of NSP1 on translation.

Brief considerations on targeting RNA with small molecules.

For more than three decades, RNA has been known to be a relevant and attractive macromolecule to target but figuring out which RNA should be targeted and how remains challenging. Recent years have seen the confluence of approaches for screening, drug optimization, and target validation that have led to the approval of a few RNA-targeting therapeutics for clinical applications. This focused perspective aims to highlight - but not exhaustively review - key factors accounting for these successes while pointing at crucial aspects worth considering for further breakthroughs.

Evaluation of the stereochemical quality of predicted RNA 3D models in the RNA-Puzzles submissions.

In silico prediction is a well-established approach to derive a general shape of an RNA molecule based on its sequence or secondary structure. This paper reports an analysis of the stereochemical quality of the RNA three-dimensional models predicted using dedicated computer programs. The stereochemistry of 1052 RNA 3D structures, including 1030 models predicted by fully automated and human-guided approaches within 22 RNA-Puzzles challenges and reference structures, is analyzed. The evaluation is based on standards of RNA stereochemistry that the Protein Data Bank requires from deposited experimental structures. Deviations from standard bond lengths and angles, planarity, or chirality are quantified. A reduction in the number of such deviations should help in the improvement of RNA 3D structure modeling approaches.

Bacterial translation machinery for deliberate mistranslation of the genetic code.

Inaccurate expression of the genetic code, also known as mistranslation, is an emerging paradigm in microbial studies. Growing evidence suggests that many microbial pathogens can deliberately mistranslate their genetic code to help invade a host or evade host immune responses. However, discovering different capacities for deliberate mistranslation remains a challenge because each group of pathogens typically employs a unique mistranslation mechanism. In this study, we address this problem by studying duplicated genes of aminoacyl-transfer RNA (tRNA) synthetases. Using bacterial prolyl-tRNA synthetase (ProRS) genes as an example, we identify an anomalous ProRS isoform, ProRSx, and a corresponding tRNA, tRNAProA, that are predominately found in plant pathogens from Streptomyces species. We then show that tRNAProA has an unusual hybrid structure that allows this tRNA to mistranslate alanine codons as proline. Finally, we provide biochemical, genetic, and mass spectrometric evidence that cells which express ProRSx and tRNAProA can translate GCU alanine codons as both alanine and proline. This dual use of alanine codons creates a hidden proteome diversity due to stochastic Ala→Pro mutations in protein sequences. Thus, we show that important plant pathogens are equipped with a tool to alter the identity of their sense codons. This finding reveals the initial example of a natural tRNA synthetase/tRNA pair for dedicated mistranslation of sense codons.

Comparative study on tertiary contacts and folding of RNase P RNAs from a psychrophilic, a mesophilic/radiation-resistant, and a thermophilic bacterium.

In most bacterial type A RNase P RNAs (P RNAs), two major loop-helix tertiary contacts (L8-P4 and L18-P8) help to orient the two independently folding S- and C-domains for concerted recognition of precursor tRNA substrates. Here, we analyze the effects of mutations in these tertiary contacts in P RNAs from three different species: (i) the psychrophilic bacterium Pseudoalteromonas translucida (Ptr), (ii) the mesophilic radiation-resistant bacterium Deinococcus radiodurans (Dra), and (iii) the thermophilic bacterium Thermus thermophilus (Tth). We show by UV melting experiments that simultaneous disruption of these two interdomain contacts has a stabilizing effect on all three P RNAs. This can be inferred from reduced RNA unfolding at lower temperatures and a more concerted unfolding at higher temperatures. Thus, when the two domains tightly interact via the tertiary contacts, one domain facilitates structural transitions in the other. P RNA mutants with disrupted interdomain contacts showed severe kinetic defects that were most pronounced upon simultaneous disruption of the L8-P4 and L18-P8 contacts. At 37°C, the mildest effects were observed for the thermostable Tth RNA. A third interdomain contact, L9-P1, makes only a minor contribution to P RNA tertiary folding. Furthermore, D. radiodurans RNase P RNA forms an additional pseudoknot structure between the P9 and P12 of its S-domain. This interaction was found to be particularly crucial for RNase P holoenzyme activity at near-physiological Mg2+ concentrations (2 mM). We further analyzed an exceptionally stable folding trap of the G,C-rich Tth P RNA.

An RNA-centric historical narrative around the Protein Data Bank.

Some of the amazing contributions brought to the scientific community by the Protein Data Bank (PDB) are described. The focus is on nucleic acid structures with a bias toward RNA. The evolution and key roles in science of the PDB and other structural databases for nucleic acids illustrate how small initial ideas can become huge and indispensable resources with the unflinching willingness of scientists to cooperate globally. The progress in the understanding of the molecular interactions driving RNA architectures followed the rapid increase in RNA structures in the PDB. That increase was consecutive to improvements in chemical synthesis and purification of RNA molecules, as well as in biophysical methods for structure determination and computer technology. The RNA modeling efforts from the early beginnings are also described together with their links to the state of structural knowledge and technological development. Structures of RNA and of its assemblies are physical objects, which, together with genomic data, allow us to integrate present-day biological functions and the historical evolution in all living species on earth.

Neocles B. Leontis (1955 - 2020).

None.